LU DECOMPOSITION OF TOTALLY NONNEGATIVE MATRICES 

K.R. GOODEARL AND T.H. LENAGAN 

Abstract. A uniqueness theorem for an LU decomposition of a totally nonnegative 
matrix is obtained. 



0. Introduction 

An m x n matrix M with entries from ffi. is said to be totally nonnegative if each of its 
minors is nonnegative. Further, such a matrix is totally positive if each of its minors is 
strictly positive. (Warning: in some texts, the terms totally positive and strictly totally 
positive are used for our terms totally nonnegative and totally positive, respectively.) 

Totally nonnegative matrices arise in many areas of mathematics and there has been 
considerable interest lately in the study of these matrices. For background information 
and historical references, there is the newly published book by Pinkus, [TB] and also two 
good survey articles [2] and [7]. 

In this paper, we are interested in the LU decomposition theory of totally nonnega- 
tive matrices. Cryer, [61 Theorem 1.1], has proved that any totally nonnegative matrix A 
has a decomposition A = LU with L totally nonnegative lower triangular and U totally 
nonnegative upper triangular. If, in addition, A is square and nonsingular then this de- 
composition is essentially unique, see, for example, [151 pages 50-55], especially Theorem 
2.10 and Proposition 2.11. However, in the singular case such decompositions need not be 
unique, as is pointed out in [51 Page 91]. 

The aim in this paper is to refine the methods of Cryer, [HIE], and Gasca and Pena, [9], 
to produce an LU decomposition for which there is a uniqueness result. 

A short word concerning the genesis of this result may be interesting to readers. In a 
series of recent papers, [101 [TTJ, [13] , a very close connection has emerged between the theory 
of totally nonnegative matrices and the theory of the torus invariant prime ideals of the 
algebra of quantum matrices. This opens up the possibility of using results and methods 
from one of these areas to produce results in the other. The existence of the results of 
this paper was suggested by the tensor product decomposition theorem for torus invariant 
prime ideals in quantum matrices obtained in an earlier paper of the present authors, [T21 
Theorem 3.5]. 
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Conventions. If a matrix is denoted by a given capital Roman letter, its entries will be 
denoted by the corresponding lower case letter, with subscripts. E.g., the entries of a 
matrix named L will be denoted Uj. 

When writing sets of row or column indices, we assume that the indices have been listed 
in strictly ascending order. 

Recall the standard partial order on index sets of the same cardinality, say I := {i\, . . . ,i s } 
and I' := {i[, . . . ,i' s }, where %\ < i 2 < ■ ■ • < i s and i[ < i 2 < ■ ■ ■ < i' s according to our 
convention above. Then: I < I' if and only if z*. < i' k for each k — 1, . . . , s. 

If A is a matrix and J, J are subsets of row indices and column indices for A then A(I, J) 
denotes the submatrix of A obtained by using the rows indexed by I and columns indexed 
by J. If | J | = | J|, the minor determined by A (J, J), that is, Det(A(I, J)), is denoted by 
[I\ J] a, or simply by [I\J] if there is no danger of confusion. By convention, [0\0}a '■= 1 
for any matrix A. 

1. LU DECOMPOSITION WITH SPECIFIED ECHELON FORMS 

We begin by giving an LU decomposition for certain rectangular matrices, in which the 
matrices L (respectively, U) have specified lower (respectively, upper) echelon forms. The 
specification of the matrices for which this decomposition holds, and the decomposition 
itself, hold over arbitrary fields, and we keep that generality for this section. In Section [2J 
we shall prove that all totally nonnegative real matrices satisfy the required hypotheses, 
and that for such matrices, the resulting factors L and U are also totally nonnegative (see 
Theorem 12. 1Q[) . 

1.1. Echelon forms. We say that a matrix U = (ity) is in upper echelon form (or row 
echelon form) if the following hold: 

(1) If the ith row of U is nonzero and is the leftmost nonzero entry in this row, then 
um = whenever both k > i and I < j; 

(2) If the ith row of U is zero then all the rows below it are zero. 

If, in addition to (1) and (2), there are no zero rows then we say that U is in strictly upper 
echelon form. 

Similar definitions are made for lower triangular matrices. Namely, a matrix L = (Zy) is 
in lower echelon form provided the transpose of L is in upper echelon form, that is: 

(1) If the jth column of L is nonzero and kj is the uppermost nonzero entry in this 
column, then 1^ = whenever k < i and I > j; 

(2) If the jth column of L is zero then all the columns to the right of it are zero. 

If, in addition to (1) and (2), there are no zero columns then we say that L is in strictly 
lower echelon form. 

In order to obtain the desired uniqueness results, we need to be more precise concerning 
the echelon shapes of matrices as above. Let r := {r±, r 2 , . . . , r t } and c := {ci, c 2 , ■ ■ ■ , c t }, 
where 1 < r% < r 2 < ■ ■ ■ < r t < m and 1 < c\ < ■ ■ ■ < Ct < n. 
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(1) We say that an m x t matrix L = is in the class C r provided that for all 
j = 1, . . .,t, we have l rj j ^ and kj = for all i < rj. Further, L G £* if also 
l rj j = 1 for all j. Note that all the matrices in L r are in strictly lower echelon form. 

(2) Similarly, we say that a t x n matrix U = (m^) is in the class U c provided that for 
all z = 1, ...,£, we have u iCi ^ and Uij = for all j < q. All such matrices are in 
strictly upper echelon form. 

1.2. Some classes of matrices 

Let r := {ri, . . . ,r t } and c := {ci,...,q} be subsets of {l,...,m} and {l,...,n}, 
respectively. An m x n matrix A is said to be in the class Ai r ,c provided that 

(1) Rank(A) = t; 

(2) For each s with s < t, the minor [n, r 2 , . . . , r s |ci, C2, . . . , c s ]a is nonzero; 

(3) [/|J]a = whenever |/| — \J\ — s < t and either / ^ {ri,...,r s } or J ^ 
{ci, . . . , c s }. 

Remark 1.3. It is easy to check that a matrix belongs to at most one class Ai r ,c- However, 
in general, a matrix need not belong to any such class - consider, for example, the matrix 



for each s < t. In particular, [r 1; . . . , r t \l, . . . , t]i ^ 0, so that Rank(L) = t. 

Suppose that . . . , i s } ^ {n, . . . , r s }. Then ik < for some k. Thus, any submatrix 
of the form L({i\, . . . ,i s }, J) is a lower triangular matrix with a zero in the kth position 
on the diagonal; and so [ii, . . . , i s \J]l = 0. Since all s-element index sets J C {1, . . . , t} 
satisfy J > {1, . . . , s}, we thus see that L e M. r ,[i,t]i where [1, t] := {1, . . . , t}. 

Similarly, any U G U c belongs to Ai[i t t], c , where t = \c\. 

Lemma 1.4. Suppose that L is an m x t matrix in the class C r and that U is a t x n 
matrix in the class U c . 

(i) Let s < t and let I (respectively, J) be an s-element subset of {1, ... , m} (respectively, 
{!,..., t}). Then [r 1; . . . , r s \J] L ^ i/ and only if J = {!,..., s}, and [I\J]l — if 



(ii) Let s < t and let I (respectively, J) be an s-element subset of {1, ... , t] (respectively, 
{1, . . . , n}). Then [I\ci, . . . , c s ]u ^ if and only if I = {1, . . . , s}, and [I\J]u — if 
J ^ {ci, . . . , c s }. 

(iii) A := LU is an m x n matrix in the class Ai r . c . 

Proof, (i) We already have [r 1 , . . . , r s \l, . . . , s]l ^ 0, and [L\J]l = for / ^ {r 1; . . . , r s }, 
by Remark [T73l If J = {ji, . . . , j s } and J ^ {1, . . . , s}, then some jk > k, whence r Jfc > r^. 
In this case, L({ri, . . . , r s }, J) is a lower triangular matrix whose k, k entry is zero, and so 
[n, . . .,r s \J] L = 0. 

(ii) This is proved symmetrically. 




Suppose that L G C r where r := {r%, . . . , r t }. Note that 

[ri, . . . ,r s \l, . . . , s] L = l n i ■ --Ir^ ^ 0, 



I t {r 1: ...,r s }. 
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(iii) First, Rank (A) < t, as A is the product of an m x t matrix and a t x n matrix. 
However, by the Cauchy-Binet identity (Lemma 14.31) . 

[r l} . . . } r t \c u . . . } c t ] A = [ri,...,r t |l,...,t]z,[l,...,*|ci,...,ct]i/ ^ 0, 

so Rank(A) = t. 

For any s < t, by the Cauchy-Binet identity together with (i), 

[ri,...,r s \ci,...,c s ] A = 22[r 1 ,...,r s \K) L [K\c 1 ,...,c s ) u 

K 

= [n, . . . ,r s \l, . . . ,s] L [l, . . .,s\a, . . .,c s ]u ^ 0. 

Now, suppose that we have a row index set I ^ {ri, . . . , r s }. For any s-element subset 
K of {1, . . . , t}, we have = by (i), and therefore, for any s-element subset J of 

{1, ...,n}, Lemma implies that [I\J] A = J2k j ]u = 0. Similarly, [I\J] A = 

for any /, J with \I\ — \ J\ — s < t and J {c\, . . . , c s }. Therefore A e Ai r ,c- D 

The following theorem gives an explicit LU decomposition for matrices in the classes 
Ai r ,c- Uniqueness of these decompositions will be proved once existence has been estab- 
lished. 

Theorem 1.5. Let A be an m x n matrix which belongs to the class M r ,c where r := 
{ri, ...,r t } and c := {c x , . . . ,c t }. 

Set L := (lij) and U := (uij) to be the mxt and txn matrices, respectively, with entries 
as follows: Uj := for i < rj and 

hj ■= b"i,r 2 , . . . ,r j _ 1 ,i\c 1 ,c 2 , . . ■ ,c j ] A [r 1 ,r 2 , ■ ■ ■ ,rAci,c 2 , ■ ■ -,0^ 

for i > rj, while := for j < Ci and 

:= [n,r 2 , • • • ,ri\c!,C2, ■ ■ • , Ci_i, j]A[n, r 2 , ■ ■ ■ , r^ct, c 2 , . . . , Q-i]^ 1 

for j > Ci. 

Then L belongs to the class £*, while U belongs to the class U c , and A = LU. 

Proof. It is obvious from the definitions that L e £* and U € U c ; so we need to prove that 
A = LU. The proof is by induction on min{m, n} with the cases where m = 1 or n = 1 
being trivial. In this proof, any minor [I\J] without a subscript is a minor of A; that is, 
[I\J] = [I\J] A . Minors of other matrices are given subscripts. 

Assume that m, n > 2, and suppose first that an = 0. Then either r\ > 1 or c\ > 1. It 
follows that either the first row or first column of A is zero, because A e Ai r ,c- Suppose 
that the first row of A is zero, in which case r± > 1. Let A be the (m — 1) x n matrix 
obtained from A by deleting the first row. Then A 6 M. r \c where r' := {r\ — 1, . . . , r t — 1}. 
By using the inductive hypothesis, there are matrices L, U, with entries as specified above, 
such that A = LU. Note that U = U. 
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Now, 
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and it is easy to check that 
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The case where the first column of A is zero is dealt with in a similar way. 
Next, assume that a n ^ and note that 7*1 = C\ = 1 in this case. Then, by elementary 
row operations using an as the pivot, we see that A = LA, where 
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and D = (dij) is the (m — 1) x (n — 1) matrix with entries 

dij := aj+ij+i - ai+i,i a n a ij+i = + j + ^t 1 ! 1 ]" 1 • 

Also, set £:= ([1, i + 1|1, j + 1]), so that D = [lll]- 1 ^. 

Let {ii, . . . ,i s } and {ji, . . . ,j s } be subsets of {1, . . . ,t — 1}. Then 

[ii, . . . ,i s \ji, ■ ■ -,Js}b = [l,h + 1, ■ ■ ■ ,i s + 1|1, h + l,...,j s + 1] [1 1 1]"- 1 , 
by Sylvester's identity (Lemma 14. 5p . It follows that 

[ii,---,i s \ji,---,j s ]D = [iu...,i s \ji 7 ...J s ] B [l\l]~ s 

= [l,z 1 + l,...,2 s + l|l,j 1 + l,..., Js + l][l|l]- 1 . 

From this, it follows that D belongs to the class Ai r ',c' where r' := {r2 — 1, . . . , r t — 1} and 
d := {c 2 - 1, • • • ,Ct - 1}. 

By induction, there are {m — 1) x (t — 1) and (t — 1) X (n — 1) matrices L = (lij) and 

U = (uij) such that D = LU, with = = for i < r J+ i — 1 and 

lij = fo-l,. • ■ ,rj- l,i|c 2 -l, . . . ,Cj +1 -l] D [r 2 -l, ■ ■ ■ ,r j+1 -l\c 2 -l, . . .,c j+1 -l\^ 
= [l,^2,---,^,«+l|l,C2,...,c i+ i][l|l]" 1 [l,r 2 ,...,r i+1 |l,c 2 ,...,c j+ i]" 1 [l|l] 
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for i > r-j + i — 1; while Uij = = Ui+ij+i for j < c i+ i — 1 and 

Uij = b"2~l, • • • ,r m -l|c 2 -l, • • . ,q-1, j]z)[r2-l, . . .,r— l|c 2 -l, . . . , c — l]^ 1 
= [l,r 2 ,...,r i+ i|l,c 2 , ...,c i , J 7'+l][l|l] -1 [l,r2,...,r i |l,C2, ...,Ci] _1 [l|l] 

for j > Cj+i — 1, by using (*) above. 
Now, observe that 
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From our calculations of the entries of L and U above, we see that 
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and therefore A = LU. This completes the inductive step. 
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□ 



Theorem 1.6. Suppose that A is in the class M. r ,c where r := {ri,r 2 , . . . ,r t } and c : = 
{ci, c 2 , . . . , Ct}. There are unique matrices L G C* and U EU C such that A = LU, namely 
those given in Theorem W .51 

Proof. We show that if A = (ay) is in the class A4 r c and A = LU with L = (i^) G £* 
and {7 = (w^) G W c , then the entries of L and ?7 can be uniquely specified from this 
information. The result then follows. 
Note that the equations 



(LU) ri j — l ri kUkj — lril u lj — u lj 
k=l 

specify the first row of U. Similarly, 

t 

Q'icx (7'?7)ic 1 ^ lik^kci ^il^lci 
fc=l 

for all i, which specifies the first column of L, as u\ Cl ^ 0. 
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Assume as an inductive hypothesis that the first s rows of U and the first s columns of 
L have been specified. Then 

t s 

a r s+1 j = (LU) Ts+1 j = ^ ] U- s+1 k u kj = ^ ] h- 3+1 k u kj + h- 3+1 s+l u s+l,j 
k=l k=l 

s 

^ , lr s +ik'U'kj ~\~ Us-\-lj 
k=l 

for all j, which specifies the (s + l)-st row of U because the terms in the last summation 
are already known by induction. 
Finally, 

t s 

^ica+i (-^^OjCs+i ^ ^ ^ifc^fc,c s -|-i ^ ^ ^ik^kjCs+i ^i,s+l ^s+l,c 3 +i 
fc=l fc=l 

for all i, which specifies the (s + l)-st column of L because the terms in the summation on 
the right are already known by induction and u s+ i tCs+1 ^ 0. 

This finishes the inductive step and so the result is proved. □ 

2. Modified Neville elimination 

We now restrict attention to real matrices and focus on total nonnegativity. Our aim is to 
show that any mxn totally nonnegative matrix A lies in one of the classes M. r c , and that 
the matrices L and U in the decomposition A = LU of Theorem ll.5l are totally nonnegative. 
While it is possible to prove directly that these matrices are totally nonnegative, it is 
technically much less complicated to obtain that result via a modification of the Neville 
elimination process of Gasca and Pena, [9]. 

The following is an elementary, but crucial, fact about totally nonnegative matrices. 

Lemma 2.1. Suppose that A = (ay) is a totally nonnegative mxn matrix, and that 
i < k < m and j < I < n. If = then either an = or a^j = 0. As a consequence, if 
some entry in the jth column, but below the ij entry, is nonzero then all elements in the 
ith row, but to the right of the ij entry, are also zero. 

Proof. Note that < [ik\jl] = a^aki — auakj — —an a kj', so that aua^j < 0. As au,a,kj > 
this gives the desired conclusion. □ 

First, we give an informal description of the elimination process that we will use. We 
start with a totally nonnegative matrix. The aim is to use a version of the Neville elimina- 
tion procedure to produce a final matrix in echelon form with no zero rows. If a zero row 
appears at any stage in the process then we delete it (rather than moving it to the bottom 
as in ordinary Neville elimination). Otherwise, we proceed as with Neville elimination: if 
we are clearing the lower entries in a given column and want to perform a row operation 
to replace the last nonzero entry in a column by zero, then we perform a row operation 
by subtracting a suitable multiple of the row immediately above this last position. Note 
that the entry immediately above this last position will be nonzero: this is guaranteed by 
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the above lemma. In the end we produce an upper triangular matrix U in echelon form 
which contains no zero rows. Keeping track of the operations performed produces a lower 
triangular matrix L such that A = LU. We also show that each of L and U is totally 
nonnegative. 

2.2. Invariants of the elimination algorithm 

The modified Neville algorithm starts with L := I and U := A, & totally nonnegative 
matrix, and uses two moves: (i) either delete a row of zeros of U and the corresponding 
column in L, or (ii) perform a Neville elimination move. 

The first aim is to show that at all times during the modified Neville algorithm we retain 
the totally nonnegative condition for L and U and the fact that A = LU. There are two 
moves to consider. The first deletes a row of U and the corresponding column of L. Note 
that if we delete a row or column from a totally nonnegative matrix then the new matrix 
is also totally nonnegative. 

Lemma 2.3. Let B be an m x p matrix and let C be a p x n matrix. Suppose that row 
i of C is zero. Set B' to be the m x (p — 1) matrix obtained by deleting the ith column 
of B and set C to be the {p — 1) x n matrix obtained by deleting the ith row of C . Then 
B'C = BC. 

Proof. Obvious. □ 

Lemma 2.4. (i) Suppose that r := {ri, ■ ■ ■ , r t } and that L G C r . Let L' be the matrix ob- 
tained by deleting column i from L. Then V e C r > where r' := {ri,r%, . . . , rj_i, rj+i, . . . , r t }. 

(ii) Suppose that c := {ci, C2, . . . , c t } and that U G U c . Let U' be the matrix obtained by 
deleting row i from U . Then U' G U c i where c! := {ci, c 2 , . . . , c i+ i, . . . , c t }. 

Proof. Obvious. □ 

The next results consider the effect of performing a Neville elimination move; that is, 
the row operation of subtracting a suitable multiple of row s from row s + 1 on the minors 
of a matrix of the form 
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when a st and a s+ i )t are nonzero, in order to clear the entry in position (s + l,t). Note that 
the resulting matrix has the form 

\ 



* ■ ■ 


• * 
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• * 
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a s t 
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Lemma 2.5. Suppose that A = (a^) with a st ^ and = whenever i > s and j < t. 
Suppose that a s+ i }t ^ while a s+Wtt = for all w > 1. Set B = (6^) where bij = for 
i 7^ s + 1 while 6 s +ij = a s+i,j — a s+ \ja~ta S j for all j. In particular, b s+ ij = for j < t 
while b s+ i tj = [s,s + l\tj\ A a~t for j > t. 
Then 



[I\J]a 

[I\J\a ~ a s+u a- t l [I\{s + 1} U {s}\ J] A 



Proof. Obvious from the definition of B. 



when sG/ors + l^I 
when s £ I and s + 1 E I . 



□ 



We refer to the change from A to B described in this lemma as a Neville elimination 
move. The next result shows that the totally nonnegative condition is preserved under a 
Neville elimination move. This result may be well-known, but we have been unable to find 
a clear statement in the literature. 

Proposition 2.6. In the above setting, if A is totally nonnegative then so is B. 

Proof. It follows from the fact that A is totally nonnegative and the definition of B that 
each > 0. Also, for any size minor, [I\J]b = [I\J]a > whenever s E I or s + 1 £ I. 

Suppose that I > 2 and that all minors of B of size less than I x I are > 0. Let [I\ J]b 
be an / x I minor. By the above remarks, we may assume that s ^ I and that s + 1 e I. 
Consider the following cases: 

(1) s + 1 is the least entry in /; 

(2) s + 1 E I, and there exists i £ I with % < s. 

In case (1), consider first the case where there is a j G J with j < t. Then the jih 
column of B(I, J) is zero; so [I\J]b = 0. Otherwise, note that 

a st [I\J] B = [I U {s}\ J U {t}] B = [I U {s}\ J U {t}] A > 0. 

As a st > it follows that [I\ J] B > 0. 

Next, consider case (2). If [I\{s + 1} U {s}\J] A = then [I\J] B = [I\J] A > 0, by the 
previous lemma; so we may assume that 

[I\{s + 1} U {s}\ J] B = [I\{s + 1} U {s}\ J\ A ± 0. 
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Suppose that [I\{i, s + 1} U {s}|F] B = for all subsets Y of J with \Y\ = I — 1. Then 
[I\{s + 1} U J]b = 0, by Lemma [4.21 Thus, we may assume that there exists a subset 
Y of J with \Y\ = l-l and [I\{i, s + l}U{s}\Y} B > 0. Suppose that J = YU{k}. Choose 
j 6 Y. 

Apply the Laplace relation of Lemma H~l7 a) with Ji = {j} and J2 = {j, k} while 
I — {i, s, s + 1} to obtain 

[i|j]s[s, s + l\jk] B - [s\j] B [i, s + l\jk] B + [s + l\j} B [is\jk} B = 0. 
It follows that 

[s\j] B [i, s + l\jk] B = [i\j]B[s, s + l\jk] B + [s + l\j] B [is\jk] B . 

By using Muir's law of extensible minors (Lemma 14.41) . we may introduce the I — 2 row 
indices from I\{i, s + 1} and the / — 2 column indices from Y\{j} to obtain 

[I\{i,s + l}U{s}\Y] B [I\J} B = 

[I\{s + l}\Y] B [I\{i} U {s}\J] B + [I\{i}\Y\ B [I\{s + 1} U {'Mb ■ 

Now, [I\{i, s + 1} U {s}|y]s > 0, by assumption, and each of the four minors on the right 
side of this equation is > (the two of size I — 1 by the inductive hypothesis and the two of 
size I because s is in the row set of the minor). It follows that [I|J]b > 0, as required □ 

Remark 2.7. Let E(s + l,s) be the matrix with 1 in the (s + l,s) position and zero 
elsewhere. Note that, with the above notation, 

B = (I-a a+lit iQ 1 E(s + l,s))A and A = (I + a s+1 . t a^E{s + 1, s))B. 

Note also that I + a s+1 ^ t a~^E{s + 1, s) is totally nonnegative. 

2.8. The Modified Neville Algorithm 

Let A be an m x n totally nonnegative matrix of rank t. The following algorithm outputs 
an LU decomposition of A, which, as we shall see, coincides with the one given in Theorem 

Input Set L := I, the identity m x m matrix, and U := A. Note that A = LU. 



Algorithm 

Step 1 If is in strictly upper echelon form then stop and output L and U. Otherwise, if 
there is a row of U consisting entirely of zeros, go to Step 2 and if not, then go to Step 3. 

Step 2 Suppose that L is of size mxw and U of size w xn, and that some row of U is zero. 
Choose i as large as possible so that the iih row of U is zero. Delete row i from U and 
column % from L to obtain new matrices L of size m x (w — 1) and U of size (w — 1) x n. Note 
that we still have A = LU, by Lemma [2 .3[ and that L and U are still totally nonnegative. 
Go to Step 1. 
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Step 3 Suppose that all rows of U are nonzero, but U is not in upper echelon form. By 
Lemma 12.14 the leftmost nonzero column of U must have a nonzero entry in its uppermost 
position. Set U = 

If the first column of U has two or more nonzero entries then set t = 1. Otherwise, 
set t > 1 so that the submatrix of U consisting of the first t — 1 columns is in upper 
echelon form, but that consisting of the first t columns is not. Then, in view of Lemma 
12. 1[ there is a largest integer s such that u st ,u s+ i, t ^ 0; moreover, = for i > s and 
j < t. Perform a Neville elimination move on U as in Lemma 12.51 that is, replace U by 
(/ — u s+ \jU~tE(s + l,s))U] so that in the new U we have u s +i t t = 0. At the same time, 
replace L by L(I + u s+ i^ t u~^E{s + 1, s)). Note that we still have A = LU, and that U is 
totally nonnegative by Proposition 12. 6[ while L is the product of two totally nonnegative 
matrices and so is still totally nonnegative. 

Go to Step 1. 

Theorem 2.9. The above algorithm outputs an m x t totally nonnegative matrix L £ jC* , 
for some r = {ri,r2, . . . ,rt}, and a t x n totally nonnegative matrix U £ U c , for some 
c = {ci, C2, • • • , c t }, such that A = LU. 

Proof. The algorithm outputs totally nonnegative matrices L and U, in strictly lower and 
upper echelon forms, respectively, such that A = LU. Also, note that the leading entry 
in each column of L is 1. Suppose that L £ C* and U £ U c with r = {r 1 ,r 2 , . . . ,r w } 
and c = {ci, ■ ■ ■ , c w }, As L is an m x w matrix and U is a w x n matrix, we have 
t = Rank(v4) < w. Moreover, 

[n,r 2 , . . .,r w \d,c 2 , . . -,c w ] A = [ri,r 2 , . . . ,r w |l, . . . ,w) L [l, . . .,io|ci,c 2 , . . .,c w } v ^ 0, 

by using the Cauchy-Binet identity; so t > w. Hence, w = t, as required. □ 

The above theorem, combined with the results of Section (TJ yields the main result of the 
paper: 

Theorem 2.10. Let A be an m x n totally nonnegative matrix. Then there is a unique 
pair r, c such that A £ M. r ^ c . Further, there is then a unique pair L £ £*, U £ U c such 
that A = LU. The matrices L and U are totally nonnegative. They are given explicitly in 
Theorem 11.51 

Proof. By Theorem I2.9[ there exist r, c and totally nonnegative matrices L £ £*, U £ U c 
such that A = LU, and A £ -M rjC by Lemma 11.41 As noted in Remark 11.31 r an d c are 
uniquely determined by A. The uniqueness of L and U then follows from Theorem 11.61 □ 

Theorem 11.51 and the total nonnegativity of the factors L and U are known for the case 
where A is a totally nonnegative nonsingular square matrix; see, for example, [16, Theorem 
2.10 and Proposition 2.11]. However, we have not been able to locate a prior source for 
the result just proved. 
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3. Examples 

Example 3.1. We first illustrate the modified Neville algorithm at work on the example 
considered by Cryer, [5J Page 91]. The matrix in question is 



A :-- 



Cryer exhibits two distinct LU factorisations of A into totally nonnegative factors: 



.4 



It is easy to check that A is totally nonnegative of rank one, and that A belongs to the 
class -M{ 2 },{i}. We start the algorithm with the pair {I, A}: 






A = IA 




and one can easily check that 




is (essentially) the unique decomposition of A as a product of a 3 x 1 matrix and a 1 x 3 
matrix. 

Example 3.2. A more complicated example. The algorithm reveals the class of A. 
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A = I A 



( 1 \ / 
10 




1 2 

2 4 

1 2 

\0 

/0 1 



1 \ 

2 
3 

2 / 




2 
2 



\ 



12 1 
2 




It follows that A is in A^{i i 3} i {2,4}- 

4. Appendix: Matrix identities, etc. 
For any index sets I and J, set £(I, J) := \{(i, j) <E I x J \ i > j}\. 

Lemma 4.1. (Laplace relations; see, for example, [HI pl4], [HJ eqn. (3.3.4), p26]) Let A 
be an m x n matrix, I C {1, . . . , m}, and JC {1, . . . , n}. 



(a) IfJi, J 2 C {1, . . .,n} with |Ji| + |J 2 | = then 

E (-i)"-'[/,k I ].[/^ 2 ]^(<- 1) ' , " J3,|/|JlUj2U 

hUl 2 =I { 

\I»\=\J»\ 

(b) Ifh,h C {l,...,n} wi£/i + |J 2 | = |J|, i/ien 

E [Ai^],[/ 2 i^-' i - 1) "" ; ' 1,|/lU/2|JU 

JlUj 2 = J 

1^1=1^1 



(Ji n J 2 = 0) 
(Ji n J 2 ^ 0)- 



(/i n J 2 = 0) 
(/i n/ 2 ^ 0). 



Lemma 4.2. Let A be anmxn matrix, I C {1, . . . , m}, and J C {1, . . . , n}, with \I\ — \ J\. 

(a) Fix Ji C J. If [1^ J ± ] A = for all h C I with = \ J X \, then [I\J] A = 0. 

(b) Fix h C J. // Ji] A = for all JjCJ with \ J X \ = \h\, then [I\J\ A = 0. 
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Proof. By symmetry, we need only prove (a). Set J 2 = J\Ji- There is a Laplace relation 
of the form 

[i\j] A = Yl HWMhVhU- 

hul 2 =l 

As all Ji]a = 0, by assumption, it follows that [I\ J] a = 0. □ 

Lemma 4.3. (Cauchy-Binet Identity; see, for example, [H eqn. (6), p86], [TU pl4]) Let A 
be an m x t matrix and B a t x n matrix, and let I C {1, . . . , m} and J C {1, . . . , n} be 
k- element sets with k < t. Then 

[i\j}ab = J21 i \ k U[k\J}b 

K 

where K ranges over all k- element subsets of {1, ... ,t}. 

Lemma 4.4. (Muir's Law of Extensible Minors; see, for example, [15, pl79, §187], [U 
p205]) Let F be a field and suppose that 

d 

J2c s [I s \Js][K s \L s ] = 

s=l 

is a homogeneous determinantal identity for matrices over F. Suppose that P is a set of 
row indices disjoint from each of the sets I s and Q is a set of column indices disjoint from 
each of the sets J s , with \P\ = \Q\. Then 

d 

J2cs[I s UP\J s UQ][K s UP\L s UQ] =0 

s=l 

is also a determinantal identity for matrices over F. 

Lemma 4.5. (Sylvester's Identity; see, for example, [8, p32], [3J eqn. (8), p772]) Let 
A = (aij) be an n x n matrix and let m < n. Set B = (bij) to be the (n — m) x (n — m) 
matrix where b^ := [1, . . . ,m, m + . . . ,m, m + j]. Then, 

det(B) = [1, . . . , n\l, . . . , n) A [l, . . . , m|l, . . . , m\\ m ' x . 
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